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A new application of the density matrix renormalization group (DMRG) method to a system 
composed of an interacting dot coupled to a infinite one dimensional lead is presented. This method 
enables one to study the influence of the coupling to an external lead on the thermodynamical 
properties of the dot. It is shown that this method reproduces known results for a non-interacting 
dot coupled to a lead, i.e., that for strong coupling discrete states remain in the dot. We show that 
these states are robust and do not disappear once interactions in the dot are considered. Moreover, 
due to these discrete states, Coulomb blockade affects the charging of the dot even though its 
strongly coupled to a lead. 
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INTRODUCTION 

Recently there has been much interest in the influence 
of the dot-lead couplings on the properties of a quantum 
dot. Once the dot is strongly coupled to the lead a sup- 
pression of the Coulomb blockade is predicted yj , arising 
from the fact that the number of electrons in the dot is 
no longer a good quantum number. Nevertheless, a weak 
remnant of the Coulomb blockade physics persists even 
for strong coupling [2j. 

Surprisingly, the details of the couplings to the lead 
play a very important role in determining the proper- 
ties of the dot. Konig, Gefen and Schon [3j have shown 
analytically that the tunnelling density of states of a non- 
interacting two orbital dot strongly coupled to a single 
lead contains a delta like peak. Such a feature indicates 
that a well defined localized state remains in the dot al- 
though the dot is strongly connected to the lead. Thus, if 
one measures the number of electrons in the dot as func- 
tion of the chemical potential or gate voltage, one will 
observe sharp jumps once the chemical potential crosses 
the localized state energy. Once two leads are connected 
to a two orbital dot a delta like peak in the density 
of states of the dot exists if the lead-dot coupling ma- 
trix elements are of equal signs, while no such feature 
is predicted if one of the matrix elements is of opposite 
sign. When the dot contains Nr> t orbitals connected to 
a single lead, Nr) ot — 1 delta like peaks appear in the tun- 
nelling density of states [j| . The situation for interacting 
dots has not been studied, although it was argued that 
the delta like peaks will not be suppressed by interactions 
i. 

Several numerical methods, especially exact diago- 
naliztion self consistent Hartree-Fock0 0, 0] and 
density- functional approaches^! 0], have proven them- 
selves to be a very useful tool in studying the properties 
of interacting quantum dots 12] . Unfortunately, none of 
these methods is appropriate for the study of the strong 
lead-dot coupling regime. Exact diagonalization methods 



are limited by the size of the Hilbert space. Although one 
may consider a non-interacting lead, due to the coupling 
with the interacting dot the whole system (dot + lead) 
must be treated via the full many particle states. This 
will create an exponentially large Hilbert space even for 
a moderate sized lead. Taking into account that the den- 
sity of states in the lead must be much larger than in the 
dot, its clear that exact diagonalization is not a viable op- 
tion. The treatment of larger leads using self consistent 
Hartree-Fock is possible, but at the price of using an un- 
controllable approximation which might lose important 
many body effects. Functional density methods are not 
suitable for the treatment of steep potential gradients of 
the type expected close to the dot-lead interface. Thus, 
a new numerical approach is needed. 

In this paper we shall use the density matrix renor- 
malization group method (DMRG) 

in order to treat 

a single one dimensional lead connected to a dot with 
Nrjot orbitals. We shall begin by describing the DMRG 
treatment of the lead-dot system. We shall first review 
several features of the occupation of the orbitals as func- 
tion of the chemical potential for the non-interacting case 
Finally, we shall establish that for an interacting 
dot all the Nrj ot — 1 states localized on the dot show full 
Coulomb blockade behavior in the limit of strong cou- 
pling. 

The paper is organized as follows: The second section 
contains the definition of the dot-lead model. In the third 
section we describe the details of the DMRG method de- 
veloped for the calculation. The orbital occupation as 
function of the chemical potential for the non-interacting 
and interacting dots is studied in the forth section. The 
results are discussed in the last section. 



MODEL 

The system which is discussed in this paper is a quan- 
tum dot connected to a one dimensional lead. The elec- 
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trons are treated as spinless electrons. The quantum dot 
is represented by the following Hamiltonian: 

N Dot N 
H QD = ( £i ~ V a) a \ a i + U Y1 "]<!,. (1) 

i—1 i>j 

Here a\ depicts the creation operator of an electron on 
the i-th orbital in the dot, £\ corresponds to the en- 
ergy of the orbital, Nnot is the number of orbitals, and 
U = e 2 /C, where C is the capacitance of the dot. The 
influence of an external gate coupled to the dot is taken 
into account via V g . The one dimensional lead corre- 
sponds to: 

oo 

H L ead = -t ^ C ] C 3 + 1 + h -°' ( 2 ) 

3=1 

where Cj is the creation operator of an electron on the 
j-th site of the lead, and t is the hopping matrix in the 
lead. For an infinite lead with no coupling to the dot 
this will result in a lead with a band width of At. The 
tunnelling between the dot and the lead is represented 
by: 

H DL = 2 V ia \ Cl + h.c, (3) 
»=i 

where the dot is assumed to be attached to the edge of 
the lead, and the tunnelling amplitude between the i-th 
orbital in the dot and the lead is given by Vi. 

Since the lead is assumed to be infinite one can not 
assign a fixed number of electrons in the system, as is 
customary in exact diagonalization studies of interacting 
fermionic systems. The number of electrons in our system 
is determined by the chemical potential /i. Thus, the 
Hamiltonian of the dot-lead system is governed by: 

/N Dot oc \ 

H = H Dot + H Lead + Hdl - A* 5Z a * tfl4 + X! c l c j • 

V i=l i=l J 

(4) 

Since we are interested in the number of electrons 
residing in the dot at a given chemical potential fi 
(or at a particular gate voltage V g for a given /i) at 
zero temperature, we must calculate the ground state 
eigenvector of H . In the next section we shall discuss 
the numerical methods we use to calculate this quantity. 



NUMERICAL METHOD 

As long as there are no interactions in the system (i.e., 
U = in the dot) the main obstacle for obtaining the 
ground state of H (Eq. 0J) is the infinite span of the lead. 
But since we are really interested in the number of elec- 
trons in the dot, one may expect that as long as level 



broadening of the orbitals in the dot (roughly propor- 
tional to -Kv\Vi\ 2 for weak coupling, where v is the local 
density of states in the lead) is much larger than the 
level spacing within the lead, the finite lead treatment 
of the system is accurate. Thus, replacing the infinite 
lead by a finite one with N sites makes sense as long as 
N » (At/Vi) 2 (where v = 1/4*). Taking a finite lead cou- 
pled to the dot we can exactly diagonalize H and obtain 
the ground state for any value of \x of V g . 

The situation is much more difficult when interactions 
in the dot are taken into account. Exact diagonalization 
is out of questions, since the size of the Hilbert space 
grows exponentially as function of the number of sites in 
the lead. We will use the DMRG method [H| to treat the 
system. This method is usually used for one dimensional 
interacting system and the resulting ground state is as 
accurate as the one obtained from exact diagonalization. 
In order to incorporate the dot into the DMRG method 
we propose the following procedure: First we create a 
block composed of the dot and the first site of the lead, 

/N D at \ 

H B = H Dot + H DL -// I a\a t + cjci . (5) 

and exactly diagonalize it. Then we begin an iteration 
process. An additional site (in the first iteration j = 2) 
is added to the lead, and the Hamiltonian is given by: 

Hb, = H B + c--i c j + h.c - \ic}-c r (6) 

Then a superblock Hamiltonian H BmmB n is formed by 
joining H B » to its mirror image H, b r (see Fig. [JJ, re- 
sulting in 

H BmmB R = H B . + c'jCjii + h.c + H, b r, (7) 

where j R is the index of the mirror image lead. The 
ground state ^(B • »B R ), which is a function of the 
superblock coordinates is calculated by the Lanczos 
method. Using this ground state a density matrix, 

p(Bm) — J]] ^(B • •B R )^{B • »B R ), (8) 

is formed (the trace is over the coordinates of the mirror 
reflection part of the Hamiltonian H b „ b r). The density 
matrix is diagonalized and only half of the states with 
the largest eigenvalues are retained as a truncated basis. 
A new block Hamiltonian H B is formed by projecting 
H B , onto the truncated basis. Also other relevant 
operators, such as aj, are rewritten in terms of the 
truncated basis. The new H B now replaces the one 
given in Eq. JSJ) in a new iteration cycle. These iteration 
continue until a satisfactory accuracy is obtained for the 
occupation of the orbitals in the dot. 



FIG. 1: The structure of the superblock at the forth iteration. 
Hb represents the dot and the lead. An additional site, H,, 
is added to the lead. The whole left part (Hb») is reflected 
(-H. B n)and added to the right of the block in order to form 
the superblock H b „ b r. 



RESULTS 
Non interacting dot 

We shall begin by considering the behavior of a non- 
interacting dot coupled to a one dimensional lead. We 
are interested in the behavior of the orbital occupation 
in the dot as function of the chemical potential [i or the 
gate voltage V g . As discussed in the previous section, 
we shall use numerical exact diagonalization in order to 
obtain the single electron eigenvectors (f>k and eigenvalues 
tk of the finite lead coupled to a dot system, i.e. of H 
as given in Eq. The dot's i-th orbital occupation, 

Ti-i , as function of the chemical potential, is defined in the 
following way: 



k 



\{4> k \a\ai\(j)k)\ 2 0{iJ, - efc). 



(9) 



In Fig. typical results for rii(n) as function of the cou- 
pling (a) Vi = V 2 = V and (b)Vi = 2V 2 = V for a two 
orbital dot are presented . As expected, for small values 
of V, rii shows a jump at e$, and this jump becomes less 
sharp as V increases. One might expect that for an in- 
creasing V , any signature of the original orbitals in the 
dot will be washed out. Nevertheless, quite surprising, 
that is not the case and as V increases a new jump in 
both m and n 2 appears for a value of ji in between to 
original orbital energies e\ and e 2 . The sum of the jumps 
in the occupation of both orbital is equal to one for large 
V, i.e., although the dot is strongly coupled to the lead, 
a state which contains a single electron remains localized 
in the dot. Such a behavior was indeed predicted in Ref. 
, where the tunnelling density of states for a non inter- 
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FIG. 2: The orbital occupation n\ (full lines) and n 2 (dotted 
lines) of a two orbital dot coupled to a one dimensional lead 
as function of [i for different values of the couplings Vi,V2- In 
both plots the series of curves correspond to different values 
of V\ which change between Vi = 0.05* (lowest curve at fj, = 
—2t) and Vi = t (highest curve at fj, = — 2t) in increments of 
0.05*. (a) V 2 = Vi; (b) V 2 = Vi/2. Inset: An enlargement of 
the region between /i = —0.2* and fi = 0.2*. 



acting two level system coupled to a lead was calculated. 
They show that for strong coupling the tunnelling den- 
sity of states acquires a sharp peak at an energy between 
the two levels, on the background of a very wide peak. 
This behavior is reproduced when fi is kept constant but 
the gate voltage coupled to the dot V g is changed. 

As can be seen in Fig. [21 the location in /i-space of 
the localized state depends on the ratio of V\ to V%. This 
location may be calculated using the following three state 
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Hamiltonian: 



H 



'Mi n I 



ex Vi 
e 2 V 2 
V x V 2 



(10) 



which takes into account the effect of the lead by a single 
site that is strongly coupled to both orbitals. At the 
limit of large coupling, H 3 i eve i has two eigenvalues at 
±\/Vi + V 2 , and (assuming Vx > V 2 ) a third at 



£2 + 



EX - £2 



(ii) 



The eigenvector of this state is composed exclusively from 
the two dot orbitals, and the eigenvalue i corresponds to 
the location of the jump in the orbital occupation seen in 
Fig. [21 Thus, at the limit of strong coupling, the details 
of the lead are not important in determining the energy 
of the localized state on the dot, which can be adequately 
predicted by coupling of all orbitals to a single external 
orbital 0] . This behavior is not limited to a two orbital 
dot 0] . For any number of orbitals Noot in a dot that are 
strongly coupled to the lead, the location of the Nn ot — 1 
jumps in the orbital occupation may be estimated by 
calculating the eigenvalues of an Nu ot + 1 matrix: 



/£l 
£ 2 








£2 



v 2 » 

v 3 

oJ 



(12) 



\Vi v 2 v 3 

This is shown in Fig. |3| where the derivative of the total 
orbital occupation as function of the chemical potential 
dn/d^L (where n = J2i n i) f° r a dot with ten orbitals 
is compared with the eigenvalues of the matrix given in 
Eq. (|12|) . An excellent correspondence for the position 
of the Noot — 1 peaks of dn/d/j and the position of the 
intermediate Noot — 1 eigenvalues of Eq. H12[l is obtained. 

In the following section we shall discuss the role of in- 
teractions in determining the orbital occupation of the 
dot. As previously discussed, for the interacting case we 
must use the DMRG method. It would be interesting 
nevertheless to compare the results of the DMRG ap- 
proach to exact diagonalization results. This, of course, 
is only possible for non- interacting systems. An exam- 
ple of the comparisons we performed are presented in 
Fig. 01 where the block size (i?s) was chosen to be 32 
states, resulting in the maximum matrix size H Bm , B n of 
4, 096. The iterative process continued until the occupa- 
tion of the levels in the dot changed by no more than 10 -3 
for three consecutive iterations. The correspondence be- 
tween the exact diagonalization results and the DMRG 
ones at the tails of the jumps is perfect. At the immediate 
vicinity of the jump there is a slight difference between 
the two methods. We shall return to the behavior of the 
the DMRG in the immediate vicinity of the jump in the 
next section. 
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FIG. 3: The derivative of the total orbital occupation dn/d/j. 
as function of /i for a dot with Noot = 10 described by a 
random matrix Hamiltonian and random coupling Vi between 
the dot and lead. The dashed lines correspond to the case 
where Vi — 0, the line corresponds to an average value of 
the coupling {Vi} = OAt, and the dotted line indicates the 
position of the eigenvalues of Eq. 1121 . 
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FIG. 4: The orbital occupation ni (full lines,circles) and 712 
(dotted lines, squares) of a two orbital dot coupled to a one 
dimensional lead as function of /i for different values of the 
couplings Vi = Vs. = 0.05i, 0.3i, 0.55i, 0.8t. The lines corre- 
spond to the DMRG results while the symbols correspond to 
the exact diagonalization results. 
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Interacting dot 



In this section we shall concentrate on the influence of 
electron-electron interaction in the dot on the occupation 
of levels in the dot coupled to a lead. For weak interac- 
tions we expect the usual Coulomb blockade behavior, 
i.e., the position of the jump in occupation of the j-th or- 
bital (ordered according to E\ < 62 < ■ ■ ■ < £N 0o t) should 
occur at /i = ej + (j — 1)U. Once the coupling with the 
lead becomes strong enough, it is not clear what influence 
the interactions will have on the existence of the localized 
states in the dot discussed in the previous section. 

Let us begin by discussing the role of interactions on a 
two orbital dot (ei = — l.lt, £2 = —0.9t). The orbital oc- 
cupations for weak (U — O.lt), intermediate (U = 0.2t) 
and strong (U = OAt) interactions are presented in Fig. 
|SJ The main feature of the non-interacting behavior, 
i.e., that at the limit of strong coupling between dot and 
lead a state which contains a single electron remains lo- 
calized in the dot (Fig. |2J), persists in the interacting 
case. The energy of the localized state corresponds to 
(£2 + £i)/2 + f(U), where f{U) depends on the change 
in occupation of both orbitals as the localized state is 
filled. For Vi = V 2 one can see that f(U) ~ U/2, but 
we shall analyze f(U) in detail in the next paragraph. 
Another interesting feature of interactions can be seen 
even for moderate values of Vi. The orbital occupation 
of the second orbital ri2, which is partially filled for low 
values of [i depopulate once a jump occurs in n%. This is 
an obvious response to the interaction between the two 
levels. 

Although it is not possible to directly compare the 
DMRG results for the interacting case to an exact cal- 
culation, one can nevertheless gain some insight into the 
reliability of the method by checking its sensitivity to a 
change in the block size. In Fig. [S] the DMRG results 
for a two orbital dot in which a block size of 16 states is 
compared with a block size of 32 states. Perfect corre- 
spondence between both sizes is seen at the at the tails 
of the jumps, while small deviations are seen in the im- 
mediate vicinity of the jump. This is very similar to the 
situation for the non-interacting case (Fig. 2J. Thus, the 
DMRG seems to be somewhat less accurate in the de- 
scription of the dot occupation in the immediate vicinity 
of a steep jump in occupation. 

In order to clarify the the role played by the interaction 
and coupling in determining f(U), we begin by analyz- 
ing the case for which V<z = 0. i.e., the second orbital is 
disconnected from the lead. Under such conditions the 
second orbital may acquire only two values n% = or 
n% = 1. Naively one expects that the second orbital will 
be filled once Un\{^) < [1 — 62- This is indeed the case, 
but as can be seen in Fig. Q ni(^) is not a monotonous 
function of /i, since once the second orbital is populated, 
there is a reduction in the occupation of the first orbital. 
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FIG. 5: The orbital occupation ri\ (full lines) and rii (dotted 
lines) of a two orbital interacting dot coupled to a one dimen- 
sional lead as function of /1 for (a)C7 = O.lt; (b) U = 0.2t; (c) 
U = OAt. The couplings Vi = V 2 = OMt, O.lt, . . . t. Inset: 
the total orbital occupation n = ni + ri2 as function of ji. 
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FIG. 6: The orbital occupation n\ (full lines,circles) and ri2 
(dotted lines, squares) of a two orbital dot coupled to a one di- 
mensional lead as function of \i for different values of the cou- 
plings V\ = V2 = 0.05i, 0.3t, 0.55£, 0.8t at interaction strength 
U — OAt are compared. The lines correspond to the DMRG 
results with block size of 32 states while the symbols corre- 
spond to a block size of 16. 




0.2 



1 ■ 1 ■ 1 1 1 ■ L - J — 1 1 

-1.5 -1.3 -1.1 -0.9 -0.7 -0.5 

FIG. 7: The orbital occupation ni (full lines) and ri2 (dotted 
lines) of a two orbital interacting dot (U = 0.2t) coupled to a 
one dimensional lead as function of /i for different couplings. 
In both cases Vi = 1, while the heavy lines correspond to 
V2 = and the regular lines to V2 = 1. 



DISCUSSION 



Therefore, the condition for the population of the sec- 
ond must be rewritten as Un\{^, + ) < fi + — £2, where 
fi + is the chemical potential just above the jump. Of 
course, once the coupling V\ is strong enough to com- 
pletely localize the remaining resonance, then ni(p + ) = 
rii(/i~) = 1/2 and the second orbital will populate at 
/i = £2 + U/2. The situation becomes more complicated 
once V2 7^ 0. The additional charging energy due to 
the population of the localized state is proportional to 
U(ni((i + ) — ni(/i~))(n2(/i + ) — ri2(p~)), which must be 
smaller than /1 — e in the limit of strong coupling. 

As might be expected, the situation is similar for dots 
with a higher number of orbitals. A three orbital dot is 
depicted in Fig. |S] The orbitals have energies £1 = —1.1, 
£2 = — 1,£3 = —0.9, and the interaction U = 0.2. In 
Fig. El a four orbital dot, with energies £1 = —1.2, 
£2 = —1.1, £3 = —1, £4 = —0.9, and interaction U = 0.1 
is shown. In both cases there is a transition from Njj ot 
jumps in the occupation of the dot for weak couplings 
to Njjot — 1 jumps for strong couplings. The jumps in 
the strong coupling regime are separated by U + A (A 
is the orbital spacing) as in the regular Coulomb block- 
ade case 01- The suppression of the the occupation of 
higher orbitals as lower ones become occupied is also evi- 
dent. Thus, the Coulomb blockade behavior between the 
resonances which remain localized on the dot continues 
although the dot is well connected to the external lead. 



As we have seen in the previous section, the DMRG 
method enables us to calculate the charging of a dot 
as function of the chemical potential for the whole span 
of coupling strength. Contrary to intuition, strong cou- 
pling to a single one dimensional lead does not remove 
all discrete features from the dot, as was discussed in 
Refs. S @- This conclusion remains valid also for 
a dot in which interactions between electrons are taken 
into account. These interactions manifest themselves as 
Coulomb blockade steps in the filling of the dot. 

This paper concentrated on developing the DMRG nu- 
merical method and examining the role of interaction 
at the strong coupling limit. Several important ques- 
tions remain open. The crossover from the weak coupling 
to the strong coupling regime should be quantified, and 
its dependence on the level spacing in the dot and the 
interaction studied. This is especially important since 
Ref. 0], which essentially treats the intermediate val- 
ues of coupling, sees no pronounced discrete features in 
the dot. Thus, one may speculate that an intermedi- 
ate regime of coupling, for which discrete features in the 
dot are suppressed, exists between the weak and strong 
coupling regimes in which these features are pronounced. 
Hints of this intermediate regime might be seen in Figs. 
(|8I9|) for Vi ~ 0.2, but this warrants further study. The 
influence of the dimensionality of the lead, the number 
of channels in the lead and the number of leads should 
be clarified. The DMRG method presented in this pa- 
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FIG. 8: The orbital occupation m (full lines), 712 (dotted 
lines), and 713, (dashed lines) of a three orbital interacting dot 
(U — 0.2t) coupled to a one dimensional lead as function of fi. 
The couplings Vi = V 2 = V 3 = 0.05t, O.lt, . . . 0.6t.Inset: the 
total orbital occupation n — n\ + 712 + as function of /i. 

per could, in principal, be expanded to deal with several 
one-dimensional leads, and probably also to quasi-one di- 
mensional leads. Thus, by expanding the methods intro- 
duced here, it may be possible to answer open questions 
regarding the dot-lead coupling. 
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fully acknowledged, as well as support from the Israel 
Academy of Science. 
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